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A periodic Kolmogorov type flow is implemented in a lattice gas automaton. For given aspect 
ratios of the automaton universe and within a range of Reynolds number values, the averaged flow 
evolves towards a stationary two-dimensional ABC type flow. We show the analogy between the 
streamlines of the flow in the automaton and the phase plane trajectories of a dynamical system. 
In practice flows are commonly studied by seeding the fluid with suspended particles which play 
the role of passive tracers. Since an actual flow is time-dependent and has fluctuations, the tracers 
exhibit interesting intrinsic dynamics. When tracers are implemented in the automaton and their 
trajectories are followed, we find that the tracers displacements obey a diffusion law, with "super- 
diffusion" in the direction orthogonal to the direction of the initial forcing. 

PACS numbers: 05.40.+j, 47.25.-c, 05.60,+w 



I. INTRODUCTION 

In an incompressible fluid undergoing two-dimensional flow, the equations describing the trajectories of the fluid 
elements in terms of the stream function ip take the Hamiltonian form 

x = v x = -d y ip(x, y; t) , 

y = v y = d x ifj(x,y;t) , (1) 

where v is the fluid velocity. In practice the trajectories of the fluid elements are visualized by seeding the fluid 
with suspended particles which act as passive tracers, whose velocity is given by Eqs.(^). If the flow were stationary, 
ijj ss = ip(x, y), the system would be conservative and thus integrable, and the trajectories of the tracers would coincide 
with the curves of the contour plot of the stream function. However the actual flow possesses intrinsic fluctuations, 

ip{x,y,t) = ip ss (x ) y) + 5ip(x,y;t), (2) 

which in general precludes exact analytical solution of Eqs. (jl|) . A consequence of the presence of noise is that the 
tracers eventually exhibit diffusive behavior. The problem that we address here is the evaluation of the tracer dynamics 
in the flow (|). 

II. KOLMOGOROV FLOW 

The Navier-Stokes equations for an incompressible fluid subject to an external force F read 

<9 t v+(vV)v = -p _1 VP + z/V 2 v + F , 

V-v = 0, (3) 

where P is the hydrostatic pressure, p the mass density, and v the kinematic viscosity. When the flow is two- 
dimensional and F has the form 
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F = <S>l x cos( K y) , k = (2tt/L) , (4) 

where L is the extension of the system in the y-direction, Eqs.(^) describe the periodic Kolmogorov flow (|. Alter- 
nately, with ip defined by Eqs.(Q), the Navier-Stokes equations can be recast in the form 

(d t ~ Re^V 2 )^ + J(ip, Aip) = Re-^cosiny) ; (5) 

here Re = tpo/f is the Reynolds number, where ipo is defined by tp ss = ipocos(Ky), and J(tp, Aip) is the Jacobian 

J(iP, Atp) = d x ipd y AiP - d y ipd x Aip . (6) 

By linear stability analysis, one can show that there exists a critical value of the Reynolds number where the system 
becomes unstable. This is most easily seen by considering a perturbation to the velocity field 

v x = v%* + Sv x , v y = v s y s + Sv y , (7) 

with 

= cos(k V ) ,v s y s = 0, (8) 

and evaluating the perturbation evolution in response to the imposed forcing. We write the perturbation in terms of 
its spatial Fourier components as 

5v x (x,y;t) = e tqx (p{y; t) , 5v y (x,y;t) = e iqx <p(y;i) , (9) 

with 

oo oo 

4>(y;t)= e^<p k (t), <p(y;t) = £ e**»^ fc (t). (10) 

k— — oo k— — oa 

We insert (]?]) with (^|) and ( |l0| ) into the Navier-Stokes equations which are then Laplace-Fourier transformed to obtain 
the characteristic equation and the dispersion equation; retaining only the modes k = 0, k = k, and k = —k (i.e. 
setting to zero the amplitudes of the modes with \k\ > k), we obtain from the characteristic equation the value of the 
critical Reynolds number 

Re e = V2^ { f_ + j' 1/2 , (11) 

and from the dispersion equation the modes 

Sv x (q 7 t) - e s+t , s + = -vn 2 + vq 2 e, (12) 
Sv y (q,t) — e s -' , s- = -vq 2 e 1 (13) 



with e = 1 — (Re/Re c ) . Equation (12) shows that the mode Sv y orthogonal to the direction of the forcing (along 
the x-axis) exhibits critical slowing down, i.e. s_ — * when Re —> Re c . At Re = Re c , the Kolmogorov flow becomes 
unstable, and beyond the bifurcation (Re > Re c ) goes into a 2-dimensional ABC type flow 0] with closed streamlines, 
separatrices, and infinite trajectories between separatrices, as illustrated in Fig.l. 



III. LATTICE GAS AUTOMATON DYNAMICS 

The flow to be investigated is produced by a lattice gas automaton H which we briefly describe. A lattice gas 
automaton can be viewed as a collection of particles residing in a discrete space, a regular d-dimensional lattice £, 
where they move at discrete time steps. Associated to each lattice node (with position denoted by r) there is a 
finite set of channels (labeled by Latin indices i,j,..). Each of these channels corresponds to a discrete value of the 
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velocity (cj) that a particle positioned at the specified node may have. An exclusion principle imposes that there be a 
maximum of one particle per channel. The exclusion principle is important because it allows a symbolic representation 
of the state of the system in terms of bits, and of its dynamics in terms of operations over sets of bits, which are easily 
implemented on a computer. The state of the automaton at time t is thus described by specifying the configuration 
on each and every node, i.e. the set of bits {n.j(r, i)}f =1 , r € C, for an automaton with b channels per node. The 
evolution of the automaton takes place in two stages : propagation and collision, applied sequentially at every time 
step. Particles are first moved according to their velocity: If channel i at node r is occupied, the propagation step 
displaces that particle to channel i of node r + c,Ai. This updating is done synchronously throughout the lattice. 
The second stage in the dynamics is a local collision step: As a function of the pre-collisional configuration, a new 
configuration is chosen by a prescription based on a set of (usually stochastic) rules. This step is crucial to determine 
the type of physics the automaton will exhibit at the macro- and mesoscopic levels. In particular, the collision step 
should preserve the quantities that are invariant under the dynamics of the model system. For example, it is possible 
to construct thermal automata whose rules are such that the number of particles, the momentum and the energy 
remain unchanged by the collision step Q . In general, there will be several sets of collision rules consistent with the 
invariance of the specified constants of motion under the collision step. The choice of a particular set is then dictated 
by operational convenience, or by the need to explore a particular physical regime. Here we are interested in the 
two-dimensional hydrodynamic regime and it suffices to implement a simple set of rules governing the automaton 
dynamics on a triangular lattice M. It can be shown, starting from the microscopic equations of the automaton, 
that, provided the symmetry of the lattice is sufficient, the macroscopic dynamics of the automaton is consistent 
with the Navier-Stokes equations It is one of the virtues of the lattice gas automaton that on the sole basis of 
microscopic rules in accordance with local invariance and symmetry properties, its macroscopic behavior produces 
correct hydrodynamics. 

Furthermore the lattice gas automaton exhibits two important features: (i) It possesses a large number of degrees of 
freedom; (ii) Its Boolean microscopic nature combined with stochastic micro-dynamics results in intrinsic spontaneous 
fluctuations, and it has been shown that these fluctuations capture the essentials of actual fluctuations in real fluids 
M . Therefore the lattice gas automaton can be considered as a reservoir of excitations extending over a wide range 
of frequencies and wavelengths. 



IV. TRACER DYNAMICS 



Tracers are implemented as suspended particles in the following way. The tracer, which we denote as a i-particle, 
is subjected to the cooperative effects of the fluid particles, and its dynamics results from the local dynamics of 
the automaton, that is from the combined effects of deterministic advection (due to the non-zero average velocity 
field resulting from the constraint imposed to the lattice gas) and random advection (due to the automaton intrinsic 
fluctuations). The i-particle undergoes displacements according to the average velocity of the fluid particles computed 
over a local domain of £ and over a number (/?) of time steps of the automaton ((3 At). In the absence of external 
force (Re = 0, that is ipa = in the notation of Section the dynamics of the i-particle is governed solely by 
the automaton noise and the tracers undergo random motion leading to diffusive behavior over long distances and 
long times, i.e. (8x 2 (t)) — (Sy 2 (t)) ~ 2D$t, for t At, where Do will be referred to as the coefficient of molecular 
diffusion, and At is the elementary time step of the automaton. The value of (3 can be varied to modify the effect of 
noise and thereby to tune the coefficient of molecular diffusion Dq of the i-particle. So the dynamics of the tracers 
can be represented in a space defined on a plane parallel to the £-plane and on a time scale set by a time increment 
equal to (3 At. 

We are now interested in the dynamics of the t-particles in a fluid subject to an external force. Therefore we 
impose a bias to the velocity field of the automaton by preparing its initial state by distributing the velocities of 
the fluid particles on each node of the lattice so that on the average we obtain a periodic velocity profile according 
to Eqs.(|^). With periodic boundary conditions imposed on the automaton universe, the shear triggered by the 
velocity bias produces a stationary flow in the form of a periodic array of vortices; at a moderately high value of the 
Reynolds number (Re/Re c ~ 3) the spatial structure of the streamlines becomes analogous to the topology of the 
two-dimensional ABC flow. In the absence of noise, the i-particles would follow exactly the streamlines, and their 
motion would be ballistic. However because intrinsic fluctuations are always present in the automaton (as they are in 
real fluids) , the tracers dynamics can be strongly perturbed, and it is only in the averaged automaton flow that their 
trajectories reflect the topology of the ABC flow, as illustrated in Fig. 2. 

When we consider the large-scale, long-time limit (t ^> L*/T4, where L* and V* are the vortex characteristic 
quantities), we can adopt a Fokker-Planck formulation for the i-particle distribution function (/), where the average 
is taken over the ensemble of realizations and, by averaging over a sufficiently large region of space (homogenization 
hypothesis [pf), we obtain a diffusion- type equation |1 
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dt(f) = D* a0 d a d p (f), (14) 

where D* = F(V*, L», Do) is the effective diffusion coefficient of the t-particles in the flow. On the basis of a power 
law assumption, we infer by dimensional analysis that [fL0| 



[D*\ = [i3o] p [i.] 1 - /i [K] 1 -" 1 (15) 

indicating that for < 1, and D sufficiently small (D < L*V*, which is the case in the lattice gas automaton), the 
elements D* aa (D* a p = D^Sap) should be larger than Dq, with an important quantitative difference between D* aa 
with [i > 0, and Dpp with fi <0. 

We performed lattice gas automaton simulations as described above, and from measurements of the mean-square 
displacements of the t-particles (5x 2 (t)) and (Sy 2 (t)), we find that over sufficiently long times they exhibit diffusive 
behavior, as shown in Fig. 3. Notice that in the absence of external force, the mean-square displacement of the t- 
particle (due to the sole effect of the automaton noise) is isotropic and corresponds to molecular diffusion (full line in 
Fig.3). 

The diffusion coefficients 



t-+co 2t 

(W)> 

2t 



Dl = D yy = lim ^» . (16) 



can be evaluated quantitatively from the explicit expression of the flow and assuming a Gaussian form for the space 
and time dependence of the second moment of the noise, which is compatible with the spectrum of the automaton 
fluctuations. We obtain fl 



n * D o v x 

where (Sv 2 .) denotes the automaton fluctuations, and v x the amplitude of the forcing in the velocity field acting on 
the t-particle: x{t) = 5v x {t) + v x cos(Ky). In the limit v x = 0, one has = Df, = Dq, and for v x ^ 0, it follows from 
(|l7|), that D\ > Dm > D . We find quantitative agreement between these predictions and our simulation data (see 
Fig.3). 



V. CONCLUDING COMMENTS 



We have presented an automaton approach to the problem of "turbulent diffusion" in a time-dependent flow with 
non-trivial average. We have shown the analogy between the stream lines in the averaged flow of the automaton 
and the phase plane trajectories of the corresponding dynamical system. For the full flow, we obtain agreement 
between the theoretical value of the diffusion coefficient of the tracers and the corresponding value computed from 
the lattice gas simulation data. The most important results are: (i) above the critical value of the Reynolds number 
(Re ~ 3i?e c ), the tracers dynamics remains diffusive, and (ii) "super-diffusion" is observed in the direction orthogonal 
to the direction of the forcing (D^_ 3> Do). Higher Reynolds number regimes and flows with non-stationary average 
are now being investigated. 
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FIG. 1: Contour plot of the stream function ip(x,y) at Re/Re c = 2.5. The streamlines shown here are obtained from the 
lattice gas automaton described in Section |[Il| and subject to a periodic forcing according to Eqs.fjs]). The horizontal plane is 
the yC-plane (dimensions are given in lattice units) and the vertical axis gives the amplitude of the stream function tp(x, y) (in 
arbitrary units). 



FIG. 2: Tracer trajectories in the flow produced by the automaton with the stream function shown in Fig.l. The tracer 
displacements have been averaged over 1000 time steps. Lattice size: 512 x 128 nodes; average particle density: 0.2 per channel; 
Reynolds number Re = 2.5 Re c ; Mach number Ma < 0.2. 
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FIG. 3: Mean-square displacement of tracers measured as a function of time. The data are the projections of the displacements 
along the a;-axis (o) and y-axis (•), i.e. respectively along and perpendicular to the direction of the forcing (F = F x (y)l x ). The 
lower full line shows the diffusive behavior generated by the automaton fluctuations in the absence of external force. (3 — 1. 
Distances are given in lattice units and time in automaton time steps. 



